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We explore the effects of the proximity to a superconduc- 
tor on the level density of a billiard for the two extreme cases 
that the classical motion in the billiard is chaotic or integrable. 
In zero magnetic field and for a uniform phase in the super- 
conductor, a chaotic billiard has an excitation gap equal to 
the Thouless energy. In contrast, an integrable (rectangular 
or circular) billiard has a reduced density of states near the 
Fermi level, but no gap. We present numerical calculations 
for both cases in support of our analytical results. For the 
chaotic case, we calculate how the gap closes as a function of 
magnetic field or phase difference. 

PACS numbers: 74.50.+r, 05.45. +b, 73.23.Ps, 74.80.Fp 



There exists a way in quantum mechanics to distin- 
guish classically chaotic systems from integrable systems, 
by looking at correlations between energy levels . In 
a billiard with integrable dynamics, on the one hand, 
the energy levels arc uncorrelated, and the spectrum has 
Poisson statistics. In a chaotic billiard, on the other 
hand, level repulsion leads to strong correlations and to 
Wigner-Dyson statistics ||. Although the level correla- 
tions are different, the mean level density does not dis- 
tinguish between chaotic and integrable billiards. 

In a recent paper Q we have shown that the proxim- 
ity to a superconductor makes it possible to distinguish 
chaotic from integrable billiards by looking at the den- 
sity of states. In a chaotic billiard, on the one hand, we 
have found (using random-matrix theory) that the cou- 
pling to a superconductor by means of a point contact 
opens a gap in the density of states of the order of the 
Thouless energy Et = NT5/2n. Here N is the number 
of transverse modes in the point contact, T is the tunnel 
probability per mode, and 6 is half the mean level spac- 
ing of the isolated billiard. In an integrable rectangular 
billiard, on the other hand, the density of states vanishes 
linearly near the Fermi level, without a gap. We have 
argued (using the Bohr-Sommerfeld approximation) that 
the absence of an excitation gap is generic for integrable 
systems. In these Proceedings we present numerical cal- 
culations for a chaotic billiard in support of the random- 
matrix theory, and we consider the effects of a magnetic 
field perpendicular to the billiard. We also present calcu- 
lations for an integrable circular billiard, both exact and 
in the Bohr-Sommerfeld approximation. 

The system studied is shown schematically in the inset 



of Fig. |l|. A billiard consisting of a normal metal (N) in 
a perpendicular magnetic field B is connected to two su- 
perconductors (Si, S2) by narrow leads, each containing 
N/2 transverse modes at the Fermi energy Ep. The order 
parameter in Si and S2 has a phase difference <fi £ [0, 7r]. 
Mode n couples to a superconductor with phase <j> n = <j)/2 
for 1 < n < N/2, 4> n = -0/2 for 1 + N/2 <n<N. For 
simplicity, we assume in this paper that there is no tun- 
nel barrier in the leads (r = 1). (The generalization to 
r 7^ 1 is straightforward Q.) 

The excitation spectrum of the billiard is discrete for 
energies < E < A, where A denotes the excitation 
gap in the bulk of the superconductors. We compute the 
spectrum starting from the scattering formulation of Ref. 

The billiard with the leads but without the super- 
conductors has an N x N scattering matrix Sq(E), at an 
energy E relative to the Fermi level. The modes in the 
leads which form the basis of So are chosen such that their 
wave functions are real at the NS interface. Evanescent 
modes in the leads are disregarded. The energy depen- 
dence of So is set by the Thouless energy Et, which is 
inversely proportional to the dwell time of an electron in 
the billiard. Andreev reflection at the superconductors 
scatters electrons (at energy E > 0) into holes (at energy 
— E), with a phase increment — <p n — arccos(i?/A). (Nor- 
mal reflection at the superconductors can be neglected if 
A <c Ep.) We assume that E and Et are both <c A, so 
that we may replace arccos(.E/A) by 7r/2, while retain- 
ing the ^-dependence of Sq. The excitation spectrum is 
obtained from the determinantal equation Q 



Det[l + So{E)e i *S^(-E)e- 



0, 



(1) 



where </> is a diagonal matrix with the phases 4> n on the 
diagonal. 

The scattering matrix So can be expressed in terms of 
the Hamiltonian matrix H of the isolated billiard |J , 

S (E) = \ -2mW\E-B + iirWW ji y 1 W. (2) 

The dimension of H is M x M, and the limit M — > 00 
will be taken later on. The M X N coupling matrix W 
has elements 



W, 



1/2 

<"" ~ I ~/f~ ) >m=l,2,...M,n = l,2,...N. 
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(3) 
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The energy 6 is one half of the mean level spacing of H. 
For the chaotic billiard, we assume that the Hermitian 
matrix H has the Pandey-Mehta distribution J^JtJ 



P{H) ex 



exp 



M{l + a 2 ) 
4A 5 



M 



J2 [(Re^) 2 +a- 2 (Im^) 2 ] ] , (4) 



with A = 2MS/tt. The parameter a e [0, 1] measures the 
strength of the time-reversal symmetry breaking. The 
relation between a and the magnetic flux <!> through a 
two-dimensional billiard (area A, no impurities, Fermi 
velocity vf) is fl||| 



Ma 2 = c($e/h) 2 hv F (A6 2 )- 1/2 , 



(5) 



with c a numerical coefficient of order unity. (For exam- 
Ia/tt for a circular billiard which is chaotic be- 



ple, c 

cause of diffuse boundary scattering j9|.) Time-reversal 
symmetry is effectively broken when Ma 2 ~ AT, which 
occurs for $ <C 7i/e. The effect of such weak magnetic 
fields on the bulk superconductor can be ignored. 

A key step J[o[ is to write the determinantal equation 
(|J) as an eigenvalue equation for an effective Hamiltonian 



Det(E-H eS )=0, 

-irWe- tc/, W T -H* 



Heff = 



(6a) 
(6b) 



The average density of states p(E) of the chaotic billiard 
results from 

p{E) =-7T- 1 lmTi-g(E + iO+), g(z) = ((z-ffeff)- 1 ), 

(7) 

where (• • •) denotes an average of H with distribution 
(f|) . The matrix Green function Q inherits the 2x2 block 
structure from the effective Hamiltonian H c ff. Because 
of this block structure, it is convenient to define the 2x2 
Green function 



G 





Gl 2 \ 




\ G 21 


G22 J 





TrGu Tr^i2 
Tr^2i Tr£ 22 



(8) 



Using the diagrammatic method of Refs. |T^] and p2| , 
adapted for a matrix Green function, we have derived a 
self-consistency equation for G. To highest order in 1/M, 
this equation reads 



A 



A I 



G=—Y 



7TW. 



XG11 
4- x\G 2 i 



z 



f- XXG12 
XG22 



(9) 



where we have abbreviated = (We t ^'W T ) nn and x = 
(1 — a 2 ) /(I + a 2 ). Eq. (^) is complemented with the 
boundary condition G — ► X/z for \z\ — >• 00. 




FIG. 1. Density of states of a chaotic billiard coupled to 
two superconductors by identical ballistic point contacts, for 
four values of the magnetic flux $ through the billiard. The 
phase difference <j) between the superconductors equals 5n/6 
and in the top and bottom panel, respectively. The curves 
are computed from Eqs. ([?]) and (^). The Thouless energy 
is given by Et = NS/2n, and the critical flux <!> c is denned 
by Eqs. (|) and @. 

We take the limit M — > 00 at fixed Ma 2 and 6, and 
assume in addition that N ^> 1. Eq. (||) then simplifies 
to 



G u = [i($/$ c ) 2 G u -iz/W] x 
[G 2 2 +G 12 /cos(</>/2)], 



G 



G 2 



Gi 2 , G 12 — 1 + G 11; 
where we have defined the critical flux $ c by 
Ma 2 = iA^($/$ c ) 2 . 



(10a) 
(10b) 



(11) 

The solution of Eq. ( pl)| ) for <p — and cf> = 5n/6 is 
plotted in Fig. |l| for several values of <&/<& c - For <f> = 
and <fi = the excitation gap equals E gap = aE^, with 
E T = N5/2ir and a = 2~ 3 / 2 {Vb - l) 5 / 2 w 0.6. The 
gap decreases with increasing flux $ or phase difference 
<j>. When (f> = 0, the gap closes at the critical flux $ c ~ 
(h/e)(NS/hvF) 1 / 2 A 1 / 4 . When <f> = it, there is no gap at 
any magnetic field. For <j> between and 7T, the gap closes 
at the flux Q c (<f>) given by 



2cos(^/2) 



1 + cos(0/2) 



1/2 



(12) 



So far we have used random-matrix theory to describe 
the chaotic system. As a test, we can compute the ex- 
act quantum mechanical density of states of a specific 
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FIG. 2. Histograms: density of states of a billiard cou- 
pled to a superconductor, computed from Eq. ([j]) and av- 
eraged over a range of Fermi energies. A chaotic Sinai bil- 
liard (top inset, solid histogram) is contrasted with an in- 
tegrable circular billiard (bottom inset, dashed histogram). 
The number of propagating transverse modes at the nor- 
mal-metal — superconductor interface (dotted line in the in- 
sets) equals N = 20 in the chaotic billiard and N = 30 in the 
circular billiard. The solid curve is the prediction from ran- 
dom-matrix theory, the dashed curve is the prediction from 
the Bohr-Sommerfeld approximation. 



billiard, coupled to a superconductor. Following Doron, 
Smilansky, and Frenkel Q we study a segment of a Sinai 
billiard, drawn to scale in the top inset of Fig. ||. The 
scattering matrix Sq (E) is determined by matching wave 
functions at the dotted line separating the billiard (area 
A) from the lead (width W) . The NS interface is also cho- 
sen at the dotted line. The density of states p(E) follows 
from Eq. (Q). We average p(E) over a small variation of 
E-p, such that the number of modes N = lnt[mvFW / wh] 
in the lead does not change. The result for N = 20 is 
shown in Fig. || (solid histogram), and is seen to agree 
quite well with the prediction from random-matrix the- 
ory (solid curve). There is no adjustable parameter in 
this comparison, the mean level spacing 8 following di- 
rectly from 8 = ttH 2 jmA. 

We now turn from a chaotic to an integrable billiard. 
In Ref . p| , we computed the density of states of a rectan- 
gular billiard, coupled to a superconductor by a narrow 
lead, and found a linearly vanishing p(E) for small E. 
Here, we present data for a circular billiard (radius R) 
in support of our claim that the absence of an excita- 
tion gap is generic for integrable billiards. The circular 
billiard considered is drawn to scale in the bottom inset 
of Fig. |[ The scattering matrix Sq(E) is again deter- 



mined by matching wave functions at the dotted line, 
which also determines the location of the NS interface. 
A wedge-shaped lead (opening angle 0) is chosen in order 
not to break the rotational symmetry (which simplifies 
the calculations). The density of states is averaged over 
a range of Fermi energies at fixed N = lnt[rnvpR6 /nti]. 
The result for N = 30 is the dashed histogram in Fig. |^. 

The density of states in the integrable billiard can 
be approximated with the Bohr-Sommerfeld quantization 
rule [H, 



PBS (E) = N / dsP^YsHE 
J o „ =0 



(n + ^)ttTlv f /s) 



(13) 

Here P(s) is the classical probability that an electron en- 
tering the billiard will exit after a path length s. We 
have calculated P(s) for the circular billiard shown in 
Fig. ^|, by generating a large number of classical trajec- 
tories. The resulting density of states pbs (E) is in good 
agreement with the quantum mechanical result in Fig. |^. 

To conclude, we have calculated the density of states of 
a chaotic Sinai billiard, connected to a superconductor. 
The result is in good agreement with the prediction from 
random-matrix theory The excitation gap closes at 
a critical flux <I> C through the billiard. In order of mag- 
nitude, $ c ~ (/i/e)(r crgodic /r dwcl i) 1 / 2 , where T dwcn is the 
mean dwell time of an electron in the billiard and T ergo dic 
is the time required to explore the entire available phase 
space. The precise value of $ c depends on the shape of 
the billiard, but the dependence on the tunnel probability 
r and the phase difference <j> is universal: 

I 1/2 



$ c (r,0) = $ c (i,o) 



2cos(</>/2) 



(14) 



_cos(0/2) - l + 2/r_ 

We have shown that the excitation gap, generic for 
chaotic billiards, is absent in an integrable circular bil- 
liard. Instead, the density of states p(E) oc E for small 
E, as in the rectangular billiard considered previously [Q. 
The agreement with the semi-classical Bohr-Sommerfeld 
approximation is better for the circular than for the rect- 
angular billiard. A rigorous semi-classical theory for this 
problem remains to be developed. 
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